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Abstract 

The magnetic ac susceptibility of high-temperature superconductors is shown 
to obey some scaling relations. We try to analyse this behavior within the 
framework of a common nonlinear response function of mixed state. The de- 
rived equations for critical current and ac susceptibility (x(T)) agree with the 

scaling relations of experimental data. 
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I. INTRODUCTION 



One of the most popular means of investigating vortex dynamics in the high-temperature 
superconductors(HTS) is the measurement of the response of the vortex system to ac fields 
[1,2]. Most experiments with high-T c superconductors (HTS) deal with thin flat sample in a 
perpendicular magnetic field. The application of a time- dependent field H(t) = HQ + h ac e~ tult 
to the sample surface results in an electric-field gradient in the sample interior {h ac is the 
ac-field amplitude and oj the angular frequency). This gives rise to a shielding current, 
which in turn exerts a Lorentz force on the vortices in the sample. The measurements of 
such ac response contain much valuable information about pinning and creep of vortices and 
turn out to be useful to test models for describing the ac losses which have to be carefully 
characterized and monitored for many applications. 

The controversy over the analysis of the ac response is often noteworthy because of the 
complicate interplay of hysteretic and eddy-current losses [1,2]. Critical state models like 
Bean model and its modifications are often used to analyze the results where the hysteretic 
losses dominate. However some problems with the amplitude dependence of ac responses 
still remain unsolved. Several common features of the amplitude dependence of the in-phase 
and out-of-phase susceptibilities have been observed for different kinds of materials [1-6]: a) 
a parallel shift of the in-phase susceptibility x~T curve with increasing h ac toward lower 
temperatures is observed; b) the onset of diamagnetism and dissipation does not appear to 
depend on h ac values; c) the out-of-phase peak shifts to lower temperatures with increasing 
h ac and broadens in the low temperature side; d) the absorption peak increases slightly when 
h ac increases, In Fig.l we show the data of high J c YBCO bulk material [5] for an example. 

To describe the behavior of the vortex system on a macroscopic scale, general meth- 
ods are the Maxwell equations combined with the materials equation of superconductors 
J(E, B,T). In present work we try to analyze the low-frequency ac susceptibility within the 
framework of a common nonlinear response function generally valid for all type-II supercon- 
ductors. 
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In next section we show some widely observed scaling relations of ac susceptibility exper- 
imental data.In section III we introduce the common nonlinear response function of mixed 
state. The equations of critical current and the ac susceptibility are derived in section 
IV. A discussion of the connection between the observed scaling relations and the nonlinear 
response function is given in section V. 



II. SCALING BEHAVIOR OF SUSCEPTIBILITY 

It is interesting to note that the different measured x" (T) curves in Fig.l can be repre- 
sented by a single curve in Fig. 2 when x"CO — x" (T*) and T — T p are scaled by x'Lak~x" ' (T*) 
and T* — T p respectively, with T* the irreversibility point and x"(T p ) = Xp eak ^ ne P e& k value 
of the out-of-phase susceptibility as 

x"(T) — x"(T*) _ T-T p 

Xl eak -X"(T*) Jl{ T*-T p } {) 

A rather more amazing fact is that the experimental data of susceptibility from different 
references [2,?,?] at different frequencies up to 26MHz can also be superimposed with the 
empirical scaling relation 

X"(T)-6(T-T P ) X "(T*) = T — T p 

x; eak -0(T-T p ) X "(T*) J2[ T*-T P } {) 

as shown in Fig. 3. 

The observed x" peak position T p can also be approximately described by an empirical 
relation 

[T P (h ac )-T*] a Kh ac (3) 

as illustrated by the inset of Fig.l, where T* is the limit of T p as h ac — > 0. This power 
law shift of T p to lower values by increasing amplitude h ac can be found for various kinds 
of materials. In Fig.4, we summarized some experimental h ac ~ T p data from different 
references in literature. 
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III. NONLINEAR RESPONSE 



The equations that describe the behavior of superconductor on a macroscopic scale are 
the Maxwell equations combined with the materials equation of superconductor J(E, B, T). 
Various models in literature suggested different specific forms of the materials equation. In 
the case of ideal type-II superconductors with negligible flux pinning, the material can be 
characterized by the linear equation 

E = Pf (B,T)J (4) 

with p/ ~ p n B/B c2 , the flux-flow resistivity as estimated by Bardeen and Stephen [7]. On 
the other hand, in nonideal type II superconductors with considerable pinning, the material 
is described by a set of equations E = B x v, v = v exp[—U(J)/kT]or 

E{J) = Jp f e- u ^ B ^' kT (5) 

The activation barrier U depends on J as well as additionally depends on the temperature 
T and magnetic field B. Different types ofU(J) have been suggested to approximate the real 
barrier, for instance, the Anderson-Kim model [8] with U (J) = U c {l — J/ J co ), the logarithmic 
barrier U(J) = U c \n(J co / J) [9] and the inverse power-law with U(J) = U c [(J co / J) m — 1]. 
[10-12] 

We find, if one makes a common modification to the different model barriers U(J) as 

U(J) - U(J P = J — Ej p s ) (6) 
then the corresponding modified materials equation 

E(J) = Jp f e- u ^' kT (7) 

leads to a common normalized form as 

y = xexp[--f(l +y - x) p ] (8) 
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with x and y the normalized current density and electric field respectively. 7is a pa- 
rameter characterizing the symmetry breaking of the pinned vortices system and p is an 
exponent. 

To show the connection of the nonlinear response function Eq. (8) with the critical- state 
model U(J), we start from the expression widely used for flux creep with the logarithmic 
barrier [9], 

£(J)= P/ Jexp[-^ln(^)] (9) 

Substituting J p = J — E{J)/ pf for the current density J in the bracket on the right-hand 
side of Eq. (9), we get 

£(J)= P/ Jexp[-^ln^] (10) 
The definition of barrier implies J co > J p . Using the approximation 

oo 

In = "(! " r l T « a(l " V-'T, (V > ^ (H) 

n=l 

finally we find Eq. (10) in the form 

\^)= 1 {l + y-xf (12) 

which is the general normalized form of the materials equation Eq. (8). Here we have 

Uc _ J _E(J) 
l = a—, x = —, y= — — (13) 
kl J c0 PfJco 

In earlier works, this materials equation for type-II superconductors has also been shown 
in connection with the Anderson-Kim model and the inverse power-law U(J) [13,14]. 

The numerical factor a in the approximation Eq.(ll) should be evaluated with con- 
sidering the limitation of sample size to the realistic barrier U(J) as discussed in Refs. 
[12-14]. Considering this limitation as a cut-off of the series in Eq. (11), we have 

* = Et± l n = c + Hn c ) 

where C is the Euler constant and Nc corresponds to the realistic cut-off of the series in 
Eq. (11). Usually a is of the order 2-4. Ignoring this limitation, one gets from Eq. (10) an 
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even simpler expression 

E{J) = p f J(^ kT 

■JcO 

or 

y/x={x- y y (14) 

with a = U c /kT, though the latter can not be used to interpret the case with small 
barrier and thermally assisted flux-flow (TAFF). In Fig.5, we show the numerical solutions 
of Eq. (8) and Eq. (14) for comparison. 

Therefore, the activation barrier U(J,B,T) in Eq.(5) can be explicitly expressed as 

U(J, B, T) = U C (B, T)F[J/J c0 (B, T)} (15) 

Incorporating it into the commonly observed scaling behavior of magnetic hysteresis 
M(H) in superconductors, it can be shown that U C {B,T) and J co (B,T) in Eq.(15) must take 
the following forms [15] 

U C (B,T) = m(T)B n 

J c0 (B,T) = X(T)B m (16) 



IV. CRITICAL CURRENT AND SUSCEPTIBILITY EQUATIONS 

The nonlinear response function Eq.(8)gives current- voltage characteristic of the form 

m = M -^f± { i + - (it) 

Where i>o is a prefactor with dimension of velocity and VqB pfj as discussed in [13]. 

Defining the critical current density J c by a certain criterion of electric field E c as E(J C ) = 
E c one finds from it the expression of the critical surface 

UB.T) = UB.TM - ))i + ^^^1 (18) 
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commonly used for the engineering calculation in applied superconductivity.In the ac 
susceptibility measurements we have E c = uh ac and the irreversibility temperature T*(B) 
is defined by the condition 

U c [B,T*(B)] = kT*(B)\nA (19) 

Ec 

With this condition, the critical surface equation (18) turns into the ohmic relation of flux-flow 
regime as 

J C [B,T*(B)] = E C / Pf (20) 
For T < T*{B), critical current density can be expressed as 

J ci B,T ) = MB,T )l l- { L^^ ) i + -^- ] (21) 

In the case where sample size is much smaller than the wave length I and the ac amplitude 
h ac <C H one can neglect the variation of local current density within a period and define 
two parameters 

L P = h ac /J c , r = L p /a (22) 

with a the radius of sample, 

Substituting Eq.(16)and Eq.(21)to Eq.(22),and considering the amplitude of electric field 

induced by h ac E c = uh ac ,we find the field and temperature dependency 

h ac _ h ac x T U C (B,T*) 1/p cuh ac ! 

r ~aJ c - a Jco[B > I)[1 { T*U C (B,T) } + p f J c0 l (23) 

It has been shown by Clem [3], the in-phase and out-of-phase permeabilities of a type-II 

superconducting cylinder can be expressed as 

»' = Mr), n" = n' Q g 2 {r) (24) 

With the scale function 
5 

0i (r) = r(l - — r),0 < r < 1 
lb 

= 1 + ffl<-2 + 2-K ) '' + ( -3; + 1 " + ^ )|sm(, ' ) + ( -4 + 4-l2» Sm29+ <- 

2 

+ (-^)sm4^,r > 1 



where 9{r) = sin 1 (r 2) and 

^(r) = ^r(l-^),0<r<l 
4 1 1 

= __(1 _ _) r > 1 (26) 
311 / 2r ~ V ; 

With a maximum^" 4 * = <? 2 (?" = 1) = 0.21. 

/Xpis the dimensionless differential permeability , which increases gradually with decreasing 
temperature 

li _ [ dB^) ]H=Bo= ^ T) (2?) 

For samples with different geometry we have also the expressions similar to Eq.(25)and 
Eq.(26)for susceptibilities x'and x"but with somewhat different specific forms of gi(r) and 
# 2 (r)than Eqs.(25) and (26) [2]. 

V. DISCUSSION 

The scaling behavior of ac susceptibility mentioned in section II can be understood 
in connection with the nonlinear response function in section III. Denoting the maximum 
92 IAX ( r ) — 92(r = r p ),then from Eqs.(22)-(26)we get the equation for the out-of-phase 
susceptibility peak position T p (B)'m the form 

hac 

Tp ~ aJ c (B,T p (B)) 

-— J c ( B ^v(BW-{ T , {B)Uc{B ^ Tp{B)) ) 

__j W ^ ac 1-1 foo\ 

+ Pf J c0 (B,T p (B)) 1 1 J 

where [7 c and J c0 can be expressed as [15] 

U C (B,T) = m(T)B n oc \T*(B) - TfB n 

J c0 (B, T) = \{T)B m oc [T*(B) - T] a B m (29) 

Starting from equations (28) and (29) the widely observed scaling relations 
Eqs.(l),(2),and (3) can be naturally derived. 
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In the case of low frequency as in Fig. 1, the amplitude of electric field E c induced by the 
ac magnetic field h ac is negligibly small .Thus from Eqs.(23) and (24) one finds 

X"(T)- X "(n x " _g 2 [r(T)\ 



x; eak -x"(T*) X peak 9% 

using equations (28) and (29) we get the form 



MAX 



(30) 



x"(T)-x"(T*) MAXi-i r _ Jc(B,T p (B)) _ u max v i r _ r T*(B)-T p 
X; eak -X"(T*) ~ L52 1 92{r ~ rp J C (B,T) } - [9 > 1 92{r-r p [ Tm _ T \ } 

= [g^V^ir = r P [i - ^f-r a } (3i) 

which is just the scaling relation Eq.(l) 

In the case of radio frequency the ac losses due to flux flow is significant at high tem- 
peratures near the irreversibility line. So the last terms in the right hand sides of equations 
(21),(23)and (28) can no longer be omitted. However, the terms with the Heaviside function 
9(T — T p ) in the empirical scaling relation Eq.(2) properly substract these frequency de- 
pendent contributions from the overall critical currents and susceptibilities. Thus, again we 
see 

X(T) - 0(T - T r)x "(T') _ brX] - lgAr = (I _ T-^ r) (32) 



x; mt -e(T-T r )x"(T-) ~<" ' T--T„' 
as the experimental data from different references at different frequencies up to 26MHz are 
superimposed in Fig.3. 

The amplitude effect relation Eq.(3)can also be well understood. Since the contribution 
to critical current from pinning is dominating at temperature T = T p . Omitting the last term 
in Eq.(23), one derives from Eqs.(28)and (29)naturally the empirical relation for the peak 
position of x" 

[T p (h ac )-T*] a (xh ac (33) 



VI. SUMMARY 



We find some empirical scaling relations for the ac susceptibility of high_temperature 
superconductors. Based on the analysis of the nonlinear response function of mixed state we 
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derive the critical current and susceptibility equations which lead naturally to the observed 
scaling behavior. 
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FIGURES 

FIG. 1. x an£ i x" as functions of temperature for a YBCO sample with high J c at four values of 
hoc (A:11.2 Oe;C:2.2 Oe;D:l.l Oe) with h//c axis. Inset relation between h ac and the temperature 
at the peak[5].(/ = 337 Hz) 

FIG. 2. Scaling form of the X (T) curves in Fig.l with different h ac noted by A,B,C,and D 
respectively. 

FIG. 3. The X (T) curves of high J c YBCO bulk materia^Y&c^Ci^Oy single crystal and high 
quality YbaiCu-iO-j films at different frequencies (331Hz ~ 26Hz),c&n be superimposed when the 
susceptibility is scaled as } x „ ffi ~ x ,j\L ^ e ^ T ~ T ^}}, and the temperature is scaled as [ 1 ■ 

[X (Tp)-X (T*)0(T-T p )\ L (,J -J-p) 1 

A denote the c"(T) curves in Ref.[5] with / = 337Hz;B denote the c"(T) curves 
in Fig.5(b),6(b),8(b) of Ref.[4] with / = 26,0.1, 9MHz;C denote the c"(T) curves in 
Fig.2(a), 2(b), 6(a), 6(b) of Ref.[2]. 

FIG. 4. Relation between h ac and (1 — T p /T*) in the different experiments. T p is the temperature 
at the peak of x".D:YBCO bulk sample [5];»:Single crystal of Pri^Ceo.^CuO^-y at / = 111Hz 
and ixqH = IT [6];B:Single crystal of Pri^bCeQ^CuO^-y at / = 111Hz and hqH = 0.1T [6];T:a 
disk (diameter 1mm) of YBCO film [2];^:a rectangle (2 x 3mm 2 ) of YBCO film [2];0 

:a ring (width 

50/xm) of film YBCO film [2];A:a ring (width 25fim) of YBCO film [2]. 

FIG. 5. Numerical solutions of equation (8) (open symbols) and equation (14) (lines) for com- 
parison. 
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